Voter Model on Heterogeneous Graphs 
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We study the voter dynamics model on heterogeneous graphs. We exploit the non-conservation 
of the magnetization to characterize how consensus is reached. For a network of TV nodes with 
an arbitrary but uncorrelated degree distribution, the mean time to reach consensus Tjv scales as 
Nfii/fi2, where is the fc"^ moment of the degree distribution. For a power- law degree distribution 
Uk ~ fc-", Tn thus scales as TV for > 3, as N/\nN for i/ = 3, as Af^^'^"^)/*"-!) for 2 < < 3, 
as (InTV)^ for u — 2, and as 0{1) for u < 2. These results agree with simulation data for networks 
with both uncorrelated and correlated node degrees. 
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In this letter we study the voter model on het- 
erogeneous networks and show that its behavior is dra- 
matically different than that on regular lattices. Many 
recent studies of basic statistical mechanical models on 
heterogeneous graphs have begun to elucidate how the 
dispersity in node degree (the number of links attached 
to a node) affects critical behavior. A representative 
but incomplete set of examples include percolation 
the Ising model 0, 0, 0, IE > diffusion and random 
walks la [flj |lfl- lul, as well as the voter model itself 
[Illllllllll^ 

The voter model is perhaps the simplest and most com- 
pletely solved example of cooperative behavior. For these 
reasons, our analytical predictions for the voter model on 
heterogeneous networks should provide new insights into 
the role of underlying heterogeneity on dynamical coop- 
erative behavior. In the model, each node of a graph is 
endowed with two states - spin up and spin down. The 
evolution consists of: (i) picking a random voter; (ii) 
the selected voter adopts the state of a randomly-chosen 
neighbor. These steps are repeated until a finite system 
necessarily reaches consensus. 

One basic property of the voter model is the exit prob- 
ability, namely, the probability that the system ends with 
all spins up, i?_|_(po), as a function of the initial density 
of up spins pq. Because the mean magnetization (aver- 
aged over all realizations and all histories) is conserved on 
regular lattices, and because the only possible final states 
are consensus, E+{po) = po 0]- A second basic property 
is the mean time to reach consensus, Tjy. For regular 
lattices in d dimensions, it is known that T^r scales with 
the number of nodes N as TV^ in d = 1, as TV In TV in 
d = 2, and as TV in d > 2 0, 0|. For heterogeneous 
networks, we find that T/v grows as TV//f//X2, where pk 
is the fc"' moment of the degree distribution of the net- 
work [Eq. In contrast to lattice systems, the TV 
dependence of T/v is generally sublinear. 

To understand how dispersity in node degree affects 
voter model dynamics, we first examine the illustrative 
example of the complete bipartite graph. We then ex- 
tend this approach to determine the behavior of the voter 



model on networks with power-law degree distributions, 
but with no correlations between node degrees. Finally, 
we validate our theoretical results by simulations of the 
voter model on networks with power-law degree distribu- 
tions, both with and without node degree correlations. 
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FIG. f : The complete bipartite graph Ka,b- 

Consider the voter model on the complete bipartite 
graph KaM of N — a + b nodes that are partitioned into 
two subgraphs a and b (Fig. Each node in the a 
subgraph is connected to all nodes in the b subgraph, 
and vice versa. Let Na,b be the respective number of 
up spins on each subgraph. In an update event, these 
numbers change according to 
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(1) 



For dNa , the gain term accounts for flipping a down spin 
in subgraph a due to its interaction with an up spin in 
b, while the loss term accounts for flipping an up spin 
in subgraph a. The second equation accounts for the 
evolution of TVf,. Since the time increment for an event 
is proportional to l/{a + b), the subgraph densities pa = 
Na/a and pb = Nb/b obey pa,b = Pb,a~Pa,b, with solution 

Pa.bit) = \\pam - Pb,a(0)] e'^' + i[p,(0) + Pfc(O)]. (2) 

While the sum of the subgraph densities, Pa + Pb is 
conserved, the magnetization m — {apa + bpb)/(a + b) is 
not [l^ . However, the bias in the rate equations for pa 
and Pb drive the subgraph densities to the common value 
Pa,b{oo) — ■^[pa{0) + Pb{0)] and magnetization conserva- 
tion is restored as this final state is approached. It also 
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bears mentioning that the magnetization itself is con- 
served if the update rule is link-based Q| . 

Since Pa + Pb is conserved, the sum of the subgraph 
densities in the final state equals 2 with probability Ej^. 
Thus the exit probability is 



with solution 



E+ = ]-[pa{0)+Pbm. 



(3) 



When the initial spins on the two subgraphs are oppo- 
sitely oriented, there is an equal probability of ending 
with all spins up or all spins down, independent of the 
subgraph sizes. In the extreme case of the star graph 
ifa.i, with a » 1 up spins at the periphery and a single 
down spin at the center, there is only a 50% change that 
the system ends with all spins up. 

We now study the mean time until consensus 
T]s[{pa,Pb) - either all spins up or all spins down - as 
a function of N, pa, and pij. This consensus time satisfies 
the recursion formula 
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TniPaiPb) = 'P{Pa,Pb^ Pa±-,Pb)[TN{Pa±-;Pb)+5t\ 

a a 

+ P(Pa,P6 Pa, Pb±\)[TN{Pa, Pb±\) + 5t] 

b b 

+ P{pa; Pb ^ Pa, Pb),[TN{Pa, Pb) + St], (4) 

where St — l/{a + b) = 1/N is the time step for a sin- 
gle spin-flip attempt. For example, the first term (a 
shorthand for two contributions) accounts for flipping 
a down (up) spin in subgraph a so that pa — > Pa ± ^• 
The probability for flipping a down spin in subgraph a is 

P{pa,P6 ^ Pa + l,Pb) = ^ (1 -pa)pb, where ^(1-pa) 

is the probability to choose a down spin in subgraph a 
and Pb is the probability to choose an up spin in subgraph 
h. This equation is subject to the boundary conditions 

r^(o,o) = rAr(i,i) = 0. 

Expanding this recursion formula to second order, we 
find, after straightforward algebra. 



N5t = {pa - Pb){da - db)TNiPa,Pb) 



(5) 



-^(Pq +Pb- '2paPb) + ^dfj TN{pa,Pb) 

where di denotes a partial derivative with respect to pi. 
The first term on the right accounts for a convection that 
drives the system to equal subgraph magnetizations in a 
time of order one. Subsequently, diffusive fluctuations 
govern the ultimate approach to consensus (see Fig. [21 . 
We thus compute the consensus time by replacing the 
subgraph densities pa and pb by their common value p. 
In so doing, we ignore the initial transient for t ^ 0(1), 
during which the subgraph densities are unequal. We 
also transform the derivatives with respect to pa and pb 
in Eq. (0) to derivative with respect to p to yield 
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TN{p) = -^^[{l-p)H^-p) + plnp] (7) 

Notice that if a = 0(1) and h = 0{N) (star graph), the 
consensus time ^ ^'(1)1 while if both a and b are 
0{N), then T/v ^ 0{N), as on a complete graph. 

We now extend this approach to graphs with arbitrary 
degree distributions but without degree correlations, i.e., 
we treat all nodes with the same degree as equivalent . 
We define pk as the density of up spins in the subset of 
nodes of degree k. Similar to Eq. the recursion for 
the mean consensus time on a heterogeneous graph, with 
initial densities {pk}, is: 

TN{{pk}) = ^P(fc;T^±)[rjv(Pfc±<5fe)+5<] 
fe 

+ + (8) 

k 

where P{k; ^ ^ ±) is the probability that a spin down 
(up) on a node of degree k fiips in an update, Q is the 
probability of no fiip, and Jfc = 1/ (Nuk) is the change in 
Pk when a spin flip occurs at a site of degree k. Here Uk 
is the fraction of sites with degree k. 

Since the probability of choosing a node is 1/iV, the 
spin flip probability may be written as 



P(fc;T-±) 



(9) 



V-Sy=± 



where Axy is the adjacency matrix element between 
nodes x and y (A^y = 1 if a; and y are connected and 
Axy = otherwise). Thus the second sum is the prob- 
ability that a node x with degree k chooses a neighbor 
with spin up (down). Under the mean-fleld assumption 
that neighboring node degrees are uncorrelated, we write 
Axy as kxky I p.iN , where p,i = kuk is the mean node 
degree for the graph. That is, we replace A^y by the 
probability that an edge between node x of degree kx 
and node y of degree ky exists. Then the second sum in 
Eq. (O for spin up and spin down simplifies respectively 
to 
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njPj = w 



l-UJ. 



Ml 



Namely, we decompose the nodes y according to their 
degree, and we define uj as the average degree-weighted 
density of up spins. In this formulation, each spin of given 
sign flips with the same probability that is a function 
of the degree-weighted magnetization rather than of the 
global magnetization, as in the case for degree-regular 
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graphs. Since the first sum in Eq. © gives the density 
of down (up) spin in the subset of nodes with degree fc, 

we now write P(fc; ^ +) = 'n-i^Lu{l — pi^), and similarly, 

P(fc; + — > — ) = ?T.fc(l ~ ^)Pk- Finally, the probability 
that there is no change in a single spin flip attempt is 

Q(M) = 1 - Efe (P(fc; - - +) + P(fc; + ~- -))• 

These simplifications enable us to write Eq. (jS]) as 
= ^rife [Luil-pk)iTipk + Sk) ^Tiipk}))] 

k 

+ Y.nk[{l-oj)pk{T{pk-Sk)~T{{pkm (10) 

k 

Expanding this recursion to second order we obtain 

Nst . j:ip^-.)o,T - ^^^'^-^^''^ diT, (11) 

k k 

where dk denotes the partial derivative with respect to 
Pk- The convective terms on the right-hand side again 
drive the system to the state where pk is equal to the 
weighted magnetization lo for all k. 

To check this convective behavior, we followed the evo- 
lution of single realizations of the voter model on scale- 
free graphs with degree distribution Uk oc /c^^'^ and mean 
degree ni = Sgenerated according to the MoUoy-Reed 
(MR) model [23 ■ Each node is assigned a random num- 
ber of stubs k that is drawn from a specified degree distri- 
bution. Pairs of unlinked stubs are then randomly joined. 
This construction eliminates degree correlations between 
neighboring nodes. For the initial state, we assign all 
nodes with degree larger than pi as spin down and all 
remaining nodes as spin up. A plot of the spin up densi- 
ties pg and pii for nodes of degrees 6 and 11 versus the 
degree-weighted up-spin density shows that these "sub- 
graph" densities quickly approach equal values (Fig.|2I). 
Analogous behavior occurs on the bipartite graph and on 
scale-free networks with degree correlations. 

For long times, we thus drop the convective terms and 
set Pk — uj yk. Concomitantly, we transform the partial 
derivatives with respect to k to derivatives with respect 
to to by using dkW = Ukk/pi to reduce (|11|) to 
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(12) 



Since X^fc ~ M2: the second moment of the degree 
distribution, this equation can be reduced to a similar 
form to ©, with solution 



Tn{uj) = -iV— [{I - uj)Ml - uj) + uj\nuj]. (13) 

For a scale-free network j23| with degree distribution 
nk ~ k~^ , the vn}^ moment is pm ~ j'^'"""' k''^nk dk. Here 




FIG. 2: Trajectories of p6(t) (degree less than p\) and pi\(t) 
(degree greater than p\ = 8) versus uj, for one realiza- 
tion of the voter model on a network of 2 x 10^ nodes, 
with degree distribution nt ~ k~^'^ . The initial state is 
(Pfc>ni (0), Pfc<^i (0)) = (0, 1). The dotted curves are the ini- 
tial transient for t < 1, after which diffusive motion leads to 
consensus at (1, 1). 



k-max ^ TV^/^"^ 1) is the maximal degree in a finite net- 
work of N nodes; this is obtained from the extremal con- 
dition /j, k-" dk = iV-i [13. Thus the second mo- 
ment diverges at the upper limit for v < 3 while the first 
moment diverges for < 2. 

Assembling the results for the moments, the mean con- 
sensus time on a scale-free graph has the TV dependence 



N i^>3, 

N/\nN 1^ = 3, 

^(2.-4)/(.-i) 2<iy<3, (14) 

0(1) V <2. 
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The prediction Tjv ^ N / In N for v — ?> may explain the 
apparent power law T^v ~ j-gported in a previous 

simulation of the voter model on such a network |l4| . 

To test our predictions, we simulated the voter model 
on the MoUoy-Reed (MR) network |2Clt and on the grow- 
ing network with redirection (GNR)|23|. The GNR is 
built by adding nodes sequentially, where each new node 
attaches either to a randomly-selected node with proba- 
bility 1 — r or to the ancestor of this target with prob- 
ability r. We chose the out degree of each node to be 
4, and redirection was applied to each outgoing link of 
the new node. This construction gives a network with a 
power-law degree distribution rik oc k~'^ , with u = 1 + 
in the range (2, oo) as r is varied between and 1. 

Fig. El shows the N dependence of T/v for representa- 
tive values of the degree exponent i/ for both the MR 
network and the GNR. The results for the two networks 
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FIG. 3: Consensus time Tjv versus A'^ on scale-free MoUoy- 
Reed networks with degree distribution nt = k~" for v = 
2.3 (+), 2.5(x), and 2.7 (*). Also shown are corresponding 
results for the GNR (open symbols). Only every second data 
point is shown. The data are all based on 10 graph realiza- 
tions, each with 100 voter model realizations. The lines are 
the expected asymptotic slopes. 



with the same v are extremely close, suggesting that de- 
gree correlations have a small effect on voter model dy- 
namics. There is also curvature in the date that orig- 
inates from finite-A^ effects. Using the maximal degree 
fcmax ~ N^/'^""'^) in the definition of the moments ulti- 
mately lead to the exponent for T/v being modified by 
the corrections, for v between 2 and 3, 



d\nT 2v-A 



dhiN 



1 



1 



where a and b are of order 1 . 

For 1/ close to 2 or 3, the leading correction term decays 
slowly in A'^, causing a discrepancy between our numer- 
ics and the theory. For example, for v = 2.3 in Fig. O 
the numerical best-fit slope to the data decreases from 
0.53 to 0.48 as we successively eliminate the first 18 data 
points. This accords well with the theoretical prediction 
of 0.46 for the slope from Eq. lO. For ly = 2.5, the 
two correction terms both decay at the same rate and 
have opposite sign. Here we may expect the smallest 
corrections, as borne out by the data - the best-fit slope 
decreases from 0.680 to 0.671 as the first 18 data points 
are deleted, while the theoretical prediction for the slope 
is 2/3. The case v = 2.7 has the slowest-decaying cor- 
rection term and here we observe the largest deviation 
between simulation and theory - the slope remains in 
the range 0.77 - 0.79 as the first 18 points are deleted, 
while theory predicts a slope of 0.82. 

In summary, the voter model on heterogeneous net- 
works approaches consensus by a two-stage process of 
quick evolution to a opinion-homogeneous state followed 
by a diffusive evolution to final consensus. By neglecting 



node degree correlations, the consensus time T/v on scale- 
free graphs has the following dependence on the degree 
distribution exponent v. for v < 2, ^ C>{1), while for 
I' > 3, T]y ^ N. In the intermediate regime of 2 < < 3, 
Tn ^ Ar(2i'-4)/(>'-i), Generically, Tn grows sublinearly 
with N; that is, high-degree nodes greatly accelerate the 
approach consensus. Finally, the A^-dependence of T/v is 
virtually the same for networks without and with degree 
correlations. 
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